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The combination of searches performed by the CDF and DO collaborations at the Fermilab Teva- 
tron Collider for neutral Higgs bosons produced in association with 6 quarks is reported. The data, 
corresponding to 2.6 fb~^ of integrated luminosity at CDF and 5.2 fb'^ at DO, have been collected 
in final states containing three or more 6 jets. Upper limits are set on the cross section multiplied by 
the branching ratio varying between 44 pb and 0.7 pb in the Higgs boson mass range 90 to 300 GeV, 
assuming production of a narrow scalar boson. Significant enhancements to the production of Higgs 
bosons can be found in theories beyond the standard model, for example in supersymmetry. The re- 
sults are interpreted as upper limits in the parameter space of the minimal supersymmetric standard 
model in a benchmark scenario favoring this decay mode. 
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In the standard model (SM), electroweak symmetry 
breaking is achieved through the introduction of a single 
scalar doublet field and the existence of a single neu- 
tral scalar boson Ij-Bj is predicted. However, exten- 
sions of the SM exist with richer structure. Introduc- 
ing a second Higgs doublet, such as in Type II 2-Higgs 
Doublet Models (2HDM) j\, leads to muhiple scalar 
bosons and can give scenarios with enhanced couplings 
to down-type fermions. Supersymmetry is a plausible 
extension to the SM that introduces an additional sym- 
metry between fermions and bosons. The two Higgs bo- 
son doublets in the minimal supersymmetric standard 
model (MSSM) [1, Q lead to five physical Higgs bosons: 
three neutral (collectively denoted as (j)): h, H, and A; 
and two charged: and H~ . At leading order the 
MSSM is a Type II 2HDM model, and two parameters 
are sufficient to describe the Higgs sector. These are con- 
ventionally chosen as the ratio of the two Higgs doublet 
vacuum expectation values, tan/3, and, Ma, the mass of 
the pseudoscalar boson, A. The couplings to the down- 
type fermions are enhanced by a factor of tan /3 relative 
to those in the SM. Thus, the main decay mode is —J' 6. , 
with branching fractions near 90% at leading order (the 
remainder being mostly t^t^). While searches 
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for T-lepton-pair final states at hadron colliders are rel- 
atively insensitive to higher-order radiative corrections, 
due to cancellations between the production and decay 
processes, this is not the case for decays to b_ . Thus, 
information from decays to bb together with stringent 
limits from decays to r+r^ can constrain higher-order ef- 
fects and yield additional information about electroweak 
symmetry breaking, supersymmetry or other new physics 
beyond the SM with similar final states such as pair pro- 
duction of color octet scalars 
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Since an inclusive search for (f) ^ b_ \s difficult due 
to large multijet backgrounds, these searches rely on the 
case where the 4> boson is produced in association with 
one or more b quarks. This final state with at least three 
b quarks represents a powerful search channel, with the 
third jet providing additional suppression of the large 
multijet background at a hadron collider. 

MSSM Higgs boson production has been studied at the 
CERN LEP e+e- colhder excluding Mh,A < 93 GeV for 



all tan/3 values The CDF and DO [H-il 

collaborations extended such searches to higher masses 
and for large tan /3. The most stringent upper limits on 
tan /? in searches for neutral Higgs production for masses 
above the LEP bounds come from searches in final states 
with r leptons pair produced at the Large Hadron Col- 
lider (iMil. 



32 



The results reported in this Article make use of pp col- 
lisions at a center of mass energy of 1.96 TeV. The data, 
corresponding to integrated luminosities of 2.6 fb~^ and 
5.2 fb~^ , were collected during Run II at the Fermilab 
Tevatron Collider by the CDF and DO collaborations, re- 
spectively. The combination of searches for neutral Higgs 
bosons in final states with three or more b jets 18|, |25| is 
presented. 

The CDF and DO detectors are described in Ref. 
|. A brief outline of the reconstruction of the final 
states and event selection used in these searches fol- 
lows. A cone algorithm [s^ with a radius of ?<!. = 
v/(A?/)2 -I- (A(/7)2 = 0.7 (CDF) and 0.5 (DO), where y 
is the rapidity and (p the azimuthal angle, is used to re- 
construct jets from energy deposition in the calorimeters. 
Identification of jets arising from 6-quark fragmentation 
(^-tagging) uses an algorithm based on reconstructing 
secondary vertices of charged particles displaced from the 
pp interaction vertex at CDF [s^l and a neural network 
combining lifetime information from secondary vertices 
and the minimum distance of approach of charged parti- 
cle trajectories to the primary interaction at DO |38| . 

Events are selected at CDF with dedicated online event 
selections (triggers) requiring at least two jets and us- 
ing ^-tagging information, while at least three jets are 
required at DO. For most of the data collected at DO, 
^-tagging information is also used in the trigger decision. 
Events are selected offline with at least three jets within 
the fiducial region. This is defined by requirements on 
the momentum transverse to the proton beam direction. 
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Pt, and detector pseudorapidity, ?7dot = — ln(tan(6'/2)), 
of jets, where 9 is the polar angle with respect to the 
proton beam direction and an origin at the center of the 
detector: three jets withp^ > 20 GeV and |?7detl < 2.0 at 
CDF and three jets with px > 20 GeV and at least two 
jets with Pt > 25 GeV, |77dct| < 2.5 at DO. Events with 
more jets are accepted at CDF but only the leading (in 
Pt) four jets are used in the analysis. Exclusive channels 
with exactly three or exactly four jets are used at DO. 
The signal sample is defined by both experiments requir- 
ing at least three 6-tagged jets. In the CDF inclusive and 
the DO 4-jet channel the two leading jets and at least one 
of the third and fourth jets are required to be b tagged. 
In addition, a large sample requiring only two 6-tagged 
jets is used by both experiments to build models of the 
multijet backgrounds. 

At CDF, the background is modeled using a sum of 
templates of the invariant mass distribution of the lead- 
ing two jets, representing contributions from different 
background modes categorized according to kinematics 
and fiavor content. The templates are constructed from 
events in the double-tagged data sample where at least 
one of the leading two jets is 6-tagged. The events are 
then weighted according to the probability for at least one 
of the jets with no 6 tag to pass the tagging requirements 
under three different assumptions as to whether it arises 
from a b quark, a c quark or a light quark (or gluon). 
This results in six different mass templates. The separa- 
tion between the different components is enhanced by us- 
ing an additional jet- flavor-sensitive discriminant, Xtags, 
that makes use of kinematic properties of the charged 
particles from secondary vertices associated with each 6- 
tagged jet. 

In the DO analysis the background model is formed by 
correcting the shape of the dijet invariant mass distribu- 
tion of a data sample with two 6-tagged jets using the 
ratio of simulated multijet samples with three 6 tags and 
two b ta gs. The simulated samples are generated using 
ALPGEN [39| with showering and hadronization carried 
out using PYTHIA (40I an d detailed simulation of the de- 
tector using GEANT 4l|. Their fiavor composition, in 
terms of the relative numbers of 6, c and light jets per 
event, is determined from a simultaneous fit to the data 
across samples with differing numbers of tagged jets, dif- 
ferent 6-tagger operating points, and in small intervals 
of the scalar sum of the transverse momenta of the jets. 
The shape correction is applied as a function of the dijet 
invariant mass and the value of a likelihood-ratio discrim- 
inant, C, designed to select signal-like events in prefer- 
ence to background-like events. Only the two possible 
jet pairings from the leading jet plus either the second 
or third leading jet are considered when forming Higgs 
candidates, and the choice that gives the highest value 
of £ > 0.65 is selected. If neither pairing in the event 
is above this threshold then the event is discarded. Sys- 
tematic uncertainties are assessed to take into account 



the imperfect modeling of the background simulation. 
However, the definition of the correction as a ratio of 
distributions significantly reduces the sensitivity of the 
final model to these imperfections. Small contributions 
to the background from top-quark-pair production are 
simulated using ALPGEN and PYTHIA. Backgrounds from 
other sources, such as Z ^ bb and single-top-quark pro- 
duction are negligible. 

The efficiency for selecting signal is estimated for both 
analyses using events generated with PYTHIA. Associated 
production of a boson and a 6 quark, gb 4>b, with 
subsequent decay bb, is used to model the signal. 
The signal cross section, experimental acceptance and the 
kinematics are corrected to next-to-leading order (NLO) 



using MCFM [42|, |43|, weighting events as a function of 



the kinematics of the leading 6-quark jet not associated 
with the Higgs decay. 

Approximately 11 500 events are selected by CDF 
with at least three 6 tags and an estimated product 
of signal efficiency and acceptance varying as a func- 
tion of between 0.18% at = 90 GeV and 0.8% 
at = 200 GeV. Approximately 15 000 and 11 000 
events with at least three 6-tagged jets are accepted 
in the DO three-jet and four-jet channels, respectively. 
The corresponding signal efficiencies for a Higgs of mass 
= 200 GeV are 1.2% and 0.6%. 

The channels are combined and exclusion limits set, 
using the modified frequentist technique 44, S], with a 
log likelihood ratio (LLR) as test statistic: 



LLR= -2 In 



p{x\Hoy 



(1) 



where p represents the probabilities that the data, X, 
are drawn from the background-only {Hq), and signal- 
plus-background [Hi) hypotheses, respectively. The like- 
lihoods are constructed from the binned two-dimensional 
distribution of the dijet invariant mass versus the xtags 
discriminant for CDF and the binned one-dimensional 
dijet invariant mass distribution for DO. Theoretical pre- 
dictions of the absolute rates for the multijet backgrounds 
have large uncertainties. Therefore, additional scale fac- 
tors are applied to the background yield in the DO anal- 
ysis and the individual templates in the CDF analysis 
and introduced into the likelihood as parameters with no 
external constraints. Systematic uncertainties are intro- 
duced as either shape or normalization variations to the 
model probability density functions. These systematics 
are governed by nuisance parameters together with Gaus- 
sian constraint terms where appropriate 46|. Sources 
of uncertainty related to a common component of the 
luminosity determination and the theoretical modeling 
of the signal production are considered fully correlated. 
All other sources of uncertainty are assumed to be un- 
correlated. The modeling of the 6-jet-tagging efficiency 
and the contamination from light-quark and gluon jets 
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fake are the dominant sources of uncertainty on the back- 
ground model. These are implemented as uncertainties 
that affect the shape of the distributions entering the like- 
lihood. Additional sub-dominant uncertainties are con- 
sidered in the DO analysis arising in the modeling of the 
trigger, jet efficiency, jet energy scale and jet resolution. 
While most of those effects have a negligible impact on 
the background model, effects that are dependent on the 
differences between b, light, and gluon jets can be signif- 
icant. Dominant experimental systematic uncertainties 
on the signal model can be attributed to luminosity (6%), 
^-tagging efficiency (11-18%) and jet energy scale (2-10% 
depending on the boson mass hypothesis). 

Limits on the product of the cross section and the 
branching ratio using the LLR test statistic are extracted. 
The limits are model-independent, apart from assuming 
a single narrow Higgs boson mass peak, dominated by ex- 
perimental resolution effects. These are summarized in 
TablelH and presented in Fig. [T] The combination gives a 
sensitivity that is better than the DO expected limit alone 
by « 25% at = 100 GeV, steadily falling to < 1% 
by = 300 GeV. Excesses of events above the SM 
background expectation are observed for = 120 and 
140 GeV with significances of 2.5 standard deviations and 
2.6 standard deviations, respectively. These are driven 
by the excesses observed in the individual contributing 
analyses of 2.8 standard deviations at = 150 GeV at 
CDF and 2.5 standard deviations at Ma, = 120 GeV at 



DO. A standard convention [47[ is used to account for the 
effect that it is more likely to find a deviation (under the 
background-only) hypothesis when several mass regions 
are probed compared with only a single hypothesis. The 
significance of the excesses in the combined analysis is 
reduced to w 2 standard deviations. 

Though these limits are the key results of this search, 
it is interesting to interpret them in terms of constraints 
on benchmark models within the MSSM. As a conse- 
quence of the enhanced couplings to b quarks at large 
tan/?, the total width of the Higgs boson increases with 
tan /?. When this width becomes comparable to the ex- 
perimental resolution of 15-20% there is an impact on 
the sensitivity of the search. When interpreting the ex- 
clusion within the MSSM, the width of the Higgs boson 
and the enhancement of the product of cross section and 
branching ratio above that of the SM are calculated us- 
ing FEYNHIGGS (48|-|53|. The width is included in the 
simulation of the signal as a function of mass and tan /3 
by convoluting a relativistic Breit-Wigner function with 
the NLO cross section [i^. Additional uncertainties for 
these model-specific limits are considered that otherwise 
cancel in the model-independent limit. For comparison 



these are derived as in previous results [18|, |22|, |25| : uncer- 
tainties on the SM signal cross section are derived from 
varying the factorization and renormalization scales by a 
factor of two and from uncertainties on the parton distri- 
bution functions. The uncertainties assessed from scale 
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TABLE I: Observed and expected upper limits at the 95% 
C.L. on the product of cross section and branching ratio 
<j{gb — > 06) X BR{<j) — )■ bb), within the acceptance for the 
highest pt b quark not arising from the Higgs decay [i^l . Ex- 
pected limits are given for the median and for ±1 and ±2 
standard deviation (s.d.) variations of the background expec- 
tation. 



Tevatron, 2.6-5.2 fb"^ 

' Observed 
Expected 
Expected ±1 s.d. 
Expected ±2 s.d. 




300 
M, [GeV] 



FIG. 1: Model independent 95% C.L. upper limits on the 
product of cross section and branching ratio for the combined 
analyses, assuming a mass degeneracy between two of the 
three neutral bosons and a Higgs boson width significantly 
smaller than the experimental resolution. The dark and hght 
shaded regions (color online) correspond to the one and two 
standard deviation bands around the median expected limit. 
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variation are taken to be 10% and the parton distribution 
functions contribute an additional 3.5-13%, depending on 
the mass hypothesis. 

The masses and coupHngs of the Higgs bosons in the 
MSSM depend, in addition to tan /? and M^, on other pa- 
rameters through radiative corrections. Limits on tan5 



as a function of Ma are derived for the m}f^'^^ scenario [54 1 
that favors the hh final state, assuming a CP-conserving 
Higgs sector [HI] and a negative value of the Higgs sector 
bilinear coupling /i. Figure [2] shows exclusion limits in the 
(tan/3. Ma) plane for this scenario. Adding a further po- 
tential theoretical uncertainty on the signal cross section 
of 20%, independent of M^, would lead to an increase of 
5% in the tan (3 limit. 




— Observed 

— Expected 
Exp. ± 1 s.d. 
Exp.± 2 s.d. 



100 



150 



200 



250 300 
Ma [GeV] 



FIG. 2: 95% C.L. lower limit in the (Ma, tan/3) plane for the 
m™"^,/! = —200 GeV, including Higgs boson width effects. 
The exclusion limit obtained from the LEP experiments is 
also shown. 

In summary, the combination of results on neutral 
Higgs boson searches in multi-6-jet events from CDF and 
DO has been presented. Upper limits are set on the prod- 
uct of the cross section and branching ratio and con- 
straints are placed in the {tan (5, Ma) plane for a partic- 
ular MSSM scenario. This combination, with more than 
half the integrated luminosity from the Tevatron still to 
be analyzed, provides the most stringent limits on neu- 
tral Higgs boson production and decay in the multi-6-jet 
mode. 
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